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Abstract 

We consider a relativistic, degenerate electron gas at zero-temperature under 
the influence of a strong, uniform, static magnetic field, neglecting any form of 
interactions. Since the density of states for the electrons changes due to the presence 
of the magnetic field (which gives rise to Landau quantization), the corresponding 
equation of state also gets modified. In order to investigate the effect of very strong 
magnetic field, we focus only on systems in which a maximum of either one, two or 
three Landau level(s) is/are occupied. This is important since, if a very large number 
of Landau levels are filled, it implies a very low magnetic field strength which yields 
back Chandrasekhar's celebrated non-magnetic results. The maximum number of 
occupied Landau levels is fixed by the correct choice of two parameters, namely 
the magnetic field strength and the maximum Fermi energy of the system. We 
study the equations of state of these one-level, two-level and three-level systems and 
compare them by taking three different maximum Fermi energies. We also find the 
effect of the strong magnetic field on the mass-radius relation of the underlying star 
composed of the gas stated above. We obtain an exciting result that, it is possible 
to have an electron degenerate static star, namely magnetized white dwarfs, with a 
mass significantly greater than the Chandrasekhar limit in the range 2.3 — 2.6Mq, 
provided it has an appropriate magnetic field strength and central density. In fact, 
recent observations of peculiar Type la supernovae - SN 2006gz, SN 2007if, SN 
2009dc, SN 2003fg - seem to suggest super-Chandrasekhar-mass white dwarfs with 
masses up to 2.4 — 2.8Mq, as their most likely progenitors. Interestingly our results 
seem to lie within the observational limits. 

Keywords: degenerate Fermi gases, stellar magnetic fields, Landau levels, equations of 
state of gases, white dwarfs 
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1 Introduction 

Neutron stars are known to have high magnetic fields as large as 10^^ G or more on 
their surfaces. Several magnetic white dwarfs have also been discovered with surface fields 
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from about 10^ G to 10^ G [T], [2], [3], [5] [6] and the physics of these objects have 
also been studied from a long time [7], [S], [H]. It is likely that stronger fields exist in the 
centers of neutron stars or even white dwarfs, the limit to which is set by the scalar virial 
theorem [10] : 

2T + W + 3U + M = 0, (1) 

where T is the total kinetic (rotational) energy, W the gravitational potential energy, 11 
arises due to the internal energy and Ai the magnetic energy. Since T and 11 are both 
positive, the maximum magnetic energy can be compared to, but can never exceed, the 
gravitational energy in an equilibrium configuration. For a star of mass M and radius R 
this gives (47ri?V3)(S^„^/87r) ~ GM^R, or B^ax ~ 2 x 10\M/Mq){R/Rq)-^ G. For 
white dwarfs this limit is 10^^ G. Ostriker and Hartwick [TTj had constructed models of 
magnetic white dwarfs with magnetic field strength B ~ 10^^ G at the center but with 
a much smaller field at the surface. Thus high interior magnetic fields in white dwarfs, 
although rarely observed in nature so far, are not completely implausible. 

It was proposed by Ginzburg [12] and Woltjer [13j that the magnetic flux 0b ~ AttBR^ 
of a star is conserved during its evolution and subsequent collapse to form a remnant 
degenerate star (flux freezing phenomenon). Thus degenerate stars of small size and large 
magnetic fields are expected to be formed from parent stars which originally could have 
quite high magnetic fields of the order ~ 10^ G [11], [15]. Thus the study of such highly 
magnetized degenerate stars will help us understand the origin and evolution of stellar 
magnetic fields. 

The mass-radius relation for (non-magnetic) white dwarfs was first determined by 
Chandrasekhar [TH]. He obtained a maximum mass for stable white dwarfs, known as the 
famous Chandrasekhar limit (~ IAAMq), such that electron degeneracy pressure is just 
adequate to counteract gravitational collapse of the star. Suh et al. [T7] obtained the 
mass-radius relation for white dwarfs with B ~ 4.4 x 10^^^^^ G. They, however, worked 
in the weak field limit (thus ignoring Landau quantization) and applied Euler-MacLaurin 
expansion to the equation of state of a fully degenerate electron gas in a strong magnetic 
field [in], in order to recover the usual equation of state in absence of a magnetic field. 
They found that both the mass and radius of magnetic white dwarfs increase compared 
to non-magnetic white dwarfs, having the same central density. 

In this paper, we consider a relativistic, degenerate electron gas at zero-temperature 
under the influence of a strong, uniform, static magnetic field. We neglect any form 
of interactions between the electrons. We study the effect of strong magnetic field on 
the equation of state of the degenerate matter and consequently obtain the mass-radius 
relation for a collapsed static star that might be composed of such matter. In order to 
highlight the effect of Landau quantization of electrons due to strong magnetic field, we 
restrict our systems to have at most one, two or three Landau level(s). 

We hypothesize the possibility of existence of purely electron degenerate stars with 
extremely high magnetic fields ~ 10^^ — 10^^ G at the center, plausibly strongly magne- 
tized white dwarfs. We also investigate the possibility of such stars having a mass greater 
than the Chandrasekhar limit, in the range 2.3 — 2.6Mq. In a simple analytical frame- 
work existence of such stars has already been reported recently [18] and its astrophysical 
implications based on numerical analysis was also discussed [IH]. Interestingly recent ob- 
servations of peculiar Type la supernovae - SN 2006gz, SN 2007if, SN 2009dc, SN 2003fg 
- seem to suggest super-Chandrasekhar-mass white dwarfs as their most likely progenitors 
[20] . [2T] . These white dwarfs are believed to have masses up to 2.4 — 2.8M0. Proposed 
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mechanism by which these white dwarfs exceed the Chandrasekhar hmit is chiefly mass 
accretion from a binary companion accompanied by differential rotation [22] • This is fun- 
damentally different from what we are proposing here. However, these observations are 
quite stimulating as not only do they support the existence of super- Chandrasekhar mass 
white dwarfs, but also our results seem to lie within the observational limits. 

The paper is organized as follows. In the next section, we first recall how the equation 
of state of a cold electron degenerate gas gets modified due to the presence of a strong 
magnetic field and then state the numerical procedure followed to obtain results. Subse- 
quently in §3 we discuss the numerical results describing the nature of the equations of 
state and the mass-radius relations. In §4 we elaborate on some of the key points of this 
work, for example, the timescale of magnetic field decay, comparison with Chandrasekhar 's 
standard results, unstable branch of the mass-radius relations, the justifications for a con- 
stant magnetic field, the non-relativistic equation for hydrostatic equilibrium, neglecting 
Coulomb interactions and the anisotropy in pressure due to strong magnetic field. Finally, 
in §5 we summarize our findings with conclusions. 

2 Equation of state for a free electron gas in a strong 
magnetic field 

2.1 Basic equations 

The energy states of a free electron in a uniform magnetic field are quantized into 
what is known as Landau orbitals, which define the motion of the electron in a plane 
perpendicular to the magnetic field. On solving the Schrodinger equation in an external, 
uniform and static magnetic field directed along the z-axis, one obtains the following 
dispersion relation [23] 

2 

= z/^e+ (2) 
where quantum number u denotes the Landau level and is given by 

v = 3 + \ + (y, (3) 

when i being the principal quantum number of the Landau level (j = 0, 1, 2,...), a = 
the spin of the electron, rrie the rest mass of the electron, h the Planck's constant and 
Pz the momentum of the electron along the z-axis which may be treated as continuous 
(the motion along the field is not quantized). The cyclotron energy is fiuJc = h{eB/mec), 
where e is the charge of the electron, c the speed of light and B the magnetic field. 
Now, the electrons can become relativistic in either of the two cases: 

(i) when the density is high enough such that the mean Fermi energy of an electron 
exceeds its rest-mass energy, 

(ii) when the cyclotron energy of the electron exceeds its rest-mass energy. 

We can define a critical magnetic field strength Be from the relation fhuJc = iTLgC^, which 
gives Be = m^c^/he = 4.414 x 10^^ G. Thus in order to study the effect of a strong 
magnetic field {B > Be) on the equation of state of a relativistic, degenerate electron gas 
we have to solve the relativistic Dirac equation. We mention here that for the present 
purpose, the magnetic field considered in the electron degenerate star originates due to 
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the flux freezing plienomenon during tfie gravitational collapse of the parent star. The 
energy eigenstates in this case turn out to be [TU] 



E.,,. = \py + ^lc\l + i^^-§-)Y''- (4) 

The main effect of the magnetic field is to modify the available density of states for the 
electrons. The number of states per unit volume in an interval Ap^ for a given Landau 
level z/ is gu{eB/h'^c)Apz, where is the degeneracy that arises due to the Landau level 
splitting, such that, = 1 for i/ = and g^ = 2 for u > 1. Therefore the electron state 
density in the absence of magnetic field 



2 

7^ 



db, (5) 



has to be replaced with 

1^^9^ dp. (6) 

V 

in case of a non-zero magnetic field. 

Now in order to calculate the electron number density rZg at zero temperature, we 
have to evaluate the integral in equation ^ from p. = to Pf{t^), which is the Fermi 
momentum of the electron for the Landau level z/, to obtain [TU] 

The Fermi energy Ep of the electrons for the Landau level u is given by 

El = pf{uYc'' + my {l + 2u^). (8) 

-De 

The upper limit z/^ of the summation in equation ([7]) is derived from the condition that 
VpiyY — 0) which implies E^ > m|;C^(l + 2z/-^) and we obtain 

or 

= (10) 

2-Dd 

where ep = Ep/nieC^, is the dimensionless Fermi energy. Bp, = B/Bc, the dimensionless 
magnetic field and epmax = Epmax/^eC^, the dimensionless maximum Fermi energy of a 
system for a given Bp, and z/^. We note that z/^ is taken be the nearest lowest integer in 
equation ( [TO|l . For example, if < z/^ < 1 for a particular value of epmax and Bp, then 
the upper limit is taken to be Um = 0. 

If we define a dimensionless Fermi momentum xplu) = pp{v) /nieC, then equations ([7]) 
and (IH]) may be written as 

2Bp """^ 



Up 



^guxp{v) (11) 
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and 

= [xf(z/)2 + 1 + 2z/fiz)]'/' (12) 

or 

XFiiy) = [el-il + 2iyBD)]'^^ (13) 

where Ae = h/rrieC, is the Compton wavelength of the electron. The matter density p can 
be written as 

p = pemnne, (14) 

where pe is the mean molecular weight per electrons and mn the mass of hydrogen atom. 
The electron energy density at zero temperature is 

xp{u) 



(2vr)2A3 ^ J \meC 



m.c- £ ,.(1 + 2.5,)^ ( (TtI^) ' (15) 



(27r)2A3 



where 



z 

i^{z) = j{l+y^Y/^dy = ^zVTT^+^ln{z+VlTz^). (16) 



Then the pressure of an electron gas in a magnetic field is given by 



dUe 



'^^,n.c^Pj,,2.B.y,,[^^^^), (17) 



where 



2.2 Procedure 



ri{z) = -zVl + z^ ln(z + a/1 + z^). 

2 2 



Referring to the qualitative discussion by Lai and Shapiro [10], let Vm also denote the 
maximum number of Landau levels occupied by a cold gas of electrons in a magnetic 
field. In this case, from equation (11 01) . Vm will be the nearest highest integer. Then if 
i^m S> 1 then the Landau energy level spacing becomes a very small fraction of the Fermi 
energy and the discrete sum over v can be replaced by an integral and we get back the 
non-magnetic results. From equation ( ITOl) we see that if the magnetic field strength is 
high (for a fixed Fermi energy), i.e., Bd ^ 1, then Vm is small and the electrons are 
restricted to the lower Landau levels only. It is in this case that the magnetic field plays 
an important role in influencing the equation of state for the relativistic degenerate gas. 

Since we are investigating the effects of high magnetic field in this work, we fix z/^, 
such that it can only take values 1, 2 or 3, which we call as one- level, two-level and three- 
level system respectively. To clarify further, by one-level system we mean where only the 
ground Landau level, z/ = 0, is occupied, two-level means where both the ground and the 
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Table 1 

Parameters for the equations of state in Figure 1. 



Epmax 


Maximum Landau level (s) 


Bd 


B in units of 10^^ G 




1 


1.5 


0.066 


2meC^ 


2 


0.75 


0.033 




3 


0.5 


0.022 




1 


199.5 


8.81 


20 meC^ 


2 


99.75 


4.40 




3 


66.5 


2.94 




1 


19999.5 


882.78 


200 rrieC^ 


2 


9999.75 


441.38 




3 


6666.5 


294.26 



first {u = 1) Landau levels are occupied and three-level means where the ground, first 
and second (z/ = 2) Landau levels are occupied. Now from equation (1T0|) we see that 
once we fix Um, we obtain a fixed Bd on supplying a desired Epmax, which corresponds 
to the maximum possible density (in a star which corresponds to its central density) for 
that Urn- Hence, in our framework, the magnetic field is in accordance with the density 
of the system. We choose Ep^ax = 2meC^, 20 mgC^ and 200 m^c^ and for each we study 
the one-level, two-level and three-level systems, giving a total of 9 cases which are listed 
in Table 1. We mention here that for a given value of Ep^ax, the value of Bp, listed here 
corresponds to a lower limit. For example, when Ep^ax = 20 mgC^, Bp, with a value of 
199.5 just results in a one-level system but, if we choose any Bp > 199.5 that would also 
lead to a one-level system. 

For each of these cases we obtain the equation of state by simultaneously solving 
equations (fTT]) . (fT^ and ( IT7|) numerically from Ep = rrieC^ to Ep = Epmax, when each 
value of Ep gives one point in the Pe — p plot. In this work we choose Pe = 2 throughout. 
Figure 1 shows the equations of state for the different cases given in Table 1, which will 
be discussed in §3. 

If we are to construct the model of a strongly magnetized star made out of electron 
degenerate matter, which is approximated to be spherical in presence of constant magnetic 
field, we require to solve the following differential equation which basically comes from 
the condition of hydrostatic equilibrium [24j 

-— -— = -A^Gp, (19) 
r-^ dr \ p ar J 

where we consider P = Pg throughout this work. See, however, the Appendix in order 
to understand the effect of deviation from spherical symmetry due to the anisotropic 
effects of magnetic field, as discussed in §4.7. Since the pressure cannot be expressed 
as an analytical function of density, unlike that of Chandrasekhar's work [16], we fit the 
equation of state with the following polytropic relation: 

P = Kp^, (20) 

with different values of the adiabatic index F in different density ranges [K being a 
dimensional constant). Thus the actual equation of state is reconstructed using multiple 



6 



polytropic equations of state. One such fit is shown in Figure 1(d) and the parameters 
K and F are stated in Table 2. The motivation behind doing such a fit is the following. 
First of all with this fitting we can determine the effect of magnetic field on the adiabatic 
index of the matter. More importantly, once we use an equation of state of the form fl20l) . 
the problem essentially reduces to solving the Lane-Emden equation which arises in the 
non-magnetic case, except that in our case, K and F also carry information about the 
magnetic field in the system. 

We briefly recall here the Lane-Emden equation, since we would be referring to some 
of its solutions in the next section. We start with equations (fT9|) and (!20|) . and write 
F = 1 + ^, where n is the polytropic index. Next, two variable transformations are made 
as follows 



P = PcO", (21) 
where pc is the central density of the star and 6' is a dimensionless variable and 

r = ae, (22) 

where ^ is another dimensionless variable and a is defined as 



l-n -I 1/2 



(23) 



4nG 

which has the dimension of length. Thus using equations f l2Up . fl2ip . 0221) and 
equation f|T9|) reduces to the famous Lane-Emden equation: 

which can be solved for a given n, subjected to the following two boundary conditions: 

= 0) = 1 (25) 

and 

If n < 5, then Q falls to zero for a finite value of ^, called as ^i, which basically denotes 
the surface of the star where the pressure goes to zero (and then density becomes zero 
too in the present context). The physical radius of the star is then given by 

R = a^i. (27) 



We note that the value of n must be such that n > — 1, so that a is real in equation (123|) 
and i? > in equation (l27|l . Each value of central density pc corresponds to a particular 
value of radius R and mass M of a star. Substituting a from equation ( 123|) in equation 
(ETD we find 



l-n r-2 



Rocpc'" =Pc' . {2i 
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The mass of the spherical star (see, however, the Appendix discussing the equations for 
an oblate spheroid) is obtained by integrating the following equation: 



^ = Anr^p. (29) 
dr 



Hence, 



R ii 
M = 4:71 j r^p dr = Ana^p^ j ^^0" d^. (30) 



Again substituting a from equation fl23|) in equation fl30|) we find 

3-n 3r-4 

M oc pc'" = ' , (31) 
and then combining equations (!28|) and (!3T!) we obtain the following mass-radius relation: 

Ro^M^=M^i. (32) 

In the present work we do not make the transformations f l^T]) and fl2^ . but directly 
solve equation ( I19p (in different density regions corresponding to the particular values of 
r or n) subjected to the same boundary conditions as in equations (125|) and (!26|) which 
are written as 

p(r = 0) = p, (33) 

and 

Hence the results in equations (128|) . (!3T]) and (132|) still remain applicable to our case. A 
plot of i? as a function of M gives the mass-radius relation for the magnetic, degenerate 
static star. 

Figure 2 shows the mass-radius relations for the one-level, two-level and three-level 
systems with Ep^ax = 20 mgC^ (results for Ep^ax = 2meC^ and 200 mgC^ also show the 
same trend). Figure 3 shows a comparison between the mass-radius relations for all the 
cases stated in Table 1. All these are discussed in detail in the next section. 



3 Numerical results 

3.1 Equations of State 

Now we come to the discussions of the results obtained. We start with the equations 
of state shown in Figure 1. Let us consider the panel (b) which shows the cases with 
Epmax = 20meC^. From Table 1 we see that the one-level system (the solid line) for 
Epmax = 20 mgC^ corresponds to a magnetic field strength Bp, = 199.5, and the two- 
level and three level systems (the dotted and dashed lines respectively) correspond to 
Bp = 99.75 and 66.5 respectively. We notice that the solid curve is free of any kink, the 
dotted curve has one kink and the dashed curve has two kinks. The kinks appear when 
there is a transition from a lower Landau level to the next and they demarcate regions 
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(a) 



(b) 




(c) (d) 



120000 I 1 1 1 1 1 1 1 1 8 




Pd Pd 



Figure 1: Equations of state in a strong magnetic field (given in Table 1) for (a) Eprnax — 
2meC^, (b) Epjnax — 20 mgC^, (c) Epmax — 200 m^c^. In all three cases the solid line, 
the dotted line and the dashed lines indicate one-level, two-level and three-level systems 
respectively. In (d) the solid line is same as the dashed line in (b), but fitted with the 
dotted line by analytical formalism (see text for details). Here Pd is the pressure in units 
of 2.668 X 10^^ erg/cc and Pd is the density in units of 2 x 10^ gm/cc. 
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of the equation of state where the pressure becomes briefly independent of density. Let 
us consider the two-level systems. The portion of the equation of state below the kink 
represents the ground Landau level and the one above the kink represents the first Landau 
level. As the Fermi energy of the electrons increases, more and more electrons occupy the 
ground Landau level and both the density and pressure of the system keep increasing. 
Once the ground level is completely filled, one observes that, on increasing the Fermi 
energy of the electrons the density increases but, the pressure remains fairly constant 
for a while, after which the pressure again starts increasing with density. It is as if, the 
increase in Fermi energy during the transition is being used by the system to move to a 
higher Landau level instead of increasing the pressure. This situation seems analogous to 
that of phase transition in matter (where the temperature remains constant with respect 
to the input heat energy during the change of phase). Similar features are also seen in 
Figures 1(a) and (c) for Epmax = '2meC^ and 200 mgC^ respectively. 

Thus looking at Figure 1, we observe that in the one- level systems (sohd hnes), all 
the electrons are in the ground Landau level and hence there is no kink. In the two- 
level systems (dotted lines), as the electrons start filling up the first Landau level, a kink 
develops in the equation of state. Finally in the three- level systems (dashed lines), there 
are two kinks — the one at the lower density indicating transition to the first Landau level 
and another at the higher density indicating transition to the second Landau level. The 
value of Epmax determines the maximum density of the system and hence the positions 
of the kinks shift accordingly in Figures 1(a), (b) and (c). 

3.2 Mass-radius relations 

Next we come to Figure 2. Here we show the mass-radius relations for the one- level, 
two-level and three level systems with Epmax = 20 rrieC^ (the explanation that follows also 
holds true for the cases with Epmax = 2meC^ and 200 mgC^). Each point in the mass- 
radius curve corresponds to a star with a particular value of central density pc which is 
supplied by us as a boundary condition {Rd and Md are the dimensionless radius and 
mass of a star respectively as defined in Figure 2 caption) . 

Figure 2(a) shows the mass-radius relation for the one-level system (Bd — 199.5). We 
see that initially as Pc increases both the mass and radius increase and then at higher 
central densities the radius becomes nearly independent of the mass (we will be explaining 
later in this section as to why such a trend is observed) . We note that the last point on this 
curve has a mass ~ 2.3 Mq and radius ~ 6.4 x 10^ cm, which corresponds to the maximum 
density point of the solid curve in Figure 1(b). This denotes the density (~ 1.16 x lO"^'^ 
gm/cc) at which the ground Landau level is completely filled. Thus a star with this pc 
and a magnetic field strength of 5 = 199.5i?c = 8.81 x 10^^ G has a mass greater than 
the Chandrasekhar limit (~ 1.44 Mq for /Xg = 2). 

Figure 2(b) shows the mass-radius relation for the two- level system (Bd — 99.75). In 
this case, a maximum mass ~ 2.3Mq is reached at a radius ~ 8.9 x 10^ cm for pc ~ 4.0 x 10^ 
gm/cc, in the same way as in Figure 2(a) when the ground Landau level is completely filled. 
After which there is a turning point in the curve, from where the mass starts decreasing. 
This turning point corresponds to the kink in the corresponding equation of state (dotted 
curve in Figure 1(b)). During the transition (when 4.0x10^ gm/cc < Pc"^ 8.0x10^ gm/cc) 
from ground to first Landau level, the radius and mass both decrease with increasing pc- 
Then there is a brief range of densities (8.0 x 10^ gm/cc < pc ^ 1-2 x 10^° gm/cc) where the 
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(a) 



(b) 



Q 
Pi 




cc 




(c) 



cc 




Figure 2: Mass-radius relations with Ep^ax — 20 mgC^ for (a) one- level system, (b) two- 
level system, (c) three-level system. Here M/j is the mass of the star in units of Mq and 
Ru is the radius of the star in units of 10^ cm (the solid, dotted and dashed lines have 
the same meaning as in Figure 1). 
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Table 2 

Parameters for the fitting function of the equation of state shown in Figure 1(d) and the 
corresponding mass-radius relation shown in Figure 2(c). 



Pd in units of 2 x lO^gm/cc 


F 


r-i 


K in CGS units 


1 — n 

R oc M~ 


- 0.096 


2.9 


0.526 


4.055 


R cx 


0.096 - 0.307 


2.4 


0.714 


1.284 


R oc 


0.307 - 1.128 


2.1 


0.909 


0.914 


R oc 


1.128 - 2.117 


0.35 


-1.538 


1.105 


unstable 


2.117 - 2.956 


4/3 


3 


0.522 


M (xR^ 


2.956 - 3.842 


2.0 


1.0 


0.246 


RocM^ 


3.842 - 4.651 


0.35 


-1.538 


2.225 


unstable 


4.651 - 6.116 


4/3 


3 


0.496 


M (xR^ 


6.116 - 7.37 


2.0 


1.0 


0.142 


i?oc M° 



radius decreases as the mass remains fairly constant and ultimately at very high densities 
(1.2 X 10^ gm/cc < Pc ^ 1-5 X 10^° gm/cc) the the radius is again nearly independent of 
the mass, as in the uppermost branch. 

Figure 2(c) shows the mass-radius relation for the three-level system {Bd = 66.5). 
Here we see two turning points, denoted by the decrease of both mass and radius, which 
correspond to the two kinks in the corresponding equation of state (dashed curve in 
Figure 1(b)). The maximum mass at the first turning point ~ 2.3M0 is reached at a radius 
~ 1.1 X 10® cm for pc ~ 2.2 x 10^ gm/cc (the first kink in the equation of state). The mass 
at the second turning point is ~ 1.2Mq which has a radius ~ 6.3 x 10'' cm for pc ~ 7.6 x 10^ 
gm/cc (the second kink in the equation of state). Just after either of the turning points 
denoted by the decrease of radius and mass both, briefly the radius decreases as the mass 
remains almost constant and finally the radius becomes nearly independent of mass. The 
maximum mass ~ 2.3M0 occurs at the central density where the ground Landau level is 
completely filled and transition to the first Landau level is about to start. We observe 
that this density follows pc (three — level) < pc (two — level) < pc (one — level) as is also 
seen from the positions of the kinks in Figure 1. 

In order to explain this behavior in further detail, we resort to the Lane-Emden rela- 
tions fl28|) . (I3T!) and (1321) . From (!28l) we note that if F > 2, then R increases with pc and if 
F = 2, then R is independent of pc- From (I3T]) we note that if F > 4/3, then M increases 
with Pc and if F = 4/3, then M is independent of pc- Finally from ( |32|) we note that if 
F > 2, then R increases with M, if F = 2, then R is independent of M and if F = 4/3, 
then M is independent of R. This is exactly what is observed in Figure 2. 

Let us look again at the mass-radius relation corresponding to Figure 2(c) (also see 
Table 2). At very low densities F is ~ 3 (which is the case for non-relativistic electrons in 
the ground Landau level; Pe{= P) oc p'^; see |25]) and the value of F keeps decreasing with 
increasing density. Up to the first turning point density pc ~ 2.2 x lO^gm/cc, both the 
radius and mass keep increasing with pc and then the radius becomes nearly independent 
of mass when F ~ 2. Then F suddenly drops to the small value ~ 0.35 which marks 
the onset of the transition from ground Landau level to first Landau level. In this region 
the pressure becomes independent of density, revealing an unstable zone in the equation 
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of state (see detailed discussion in §4.5). As the density increases further, T approaches 
the relativistic value of 4/3. In this regime we see that the radius decreases slightly as 
the mass does not change significantly, as is also true for the mass-radius relation in the 
classical non-magnetic case for F = 4/3 (see Figure 4(b)). Next F takes up a value of 2 
and the radius again becomes nearly independent of mass till it reaches the second turning 
point at density pc ~ 7.6 x 10^ gm/cc. Again during the transition from the first Landau 
level to the second, F drops to 0.35, followed by values of 4/3 and 2, which have the same 
explanations as stated above. 

We also observe that for a certain range of masses in the two-level and three-level 
systems, it is possible to have multiple values of the radius for a given value of mass. For 
instance, looking at the two-level system in Figure 2(b) in the range ~ 0.6 — 1.0 Mq, we 
observe that the same value of mass corresponds to three different values of the radius. 
Let us call them Ri, R2 and R3, such that Ri > R2 > R3 and Pc(-Ri) < Pc{R2) < PciRs)- 
To explain why we observe such a behavior, we recall that Br, oi the system is such 
that Vrn is fixed for a given value of Epmax (see Table 1 for the values). In Figure 2(b), 
Bo = 99.75 ensures that the system can have at the most two Landau levels (ground and 
first), but to what extent they will be filled depends on the Fermi energy of the electrons 
(< Epmax)- For low central densities (i.e. low Fermi energy) the electrons occupy only 
the ground Landau level. In order for the electrons to start occupying the first Landau 
level, the central density of the star must be adequately high. Thus, based on previous 
discussion for the mass-radius curves in Figure 2, it is possible that the same mass of 
a star corresponds to more than one radius, depending on the Landau level occupancy. 
Now, Ri lies on the branch of the mass-radius relation which corresponds to the stars 
in which the electrons occupy only the ground Landau level, while R^ lies on the branch 
which corresponds to the electrons occupying the first Landau level (the ground level being 
already filled). R2 lies on the unstable branch of the mass-radius relation (see discussion 
in §4.5) when the electrons are in the transition mode form the ground to the first Landau 
level. The three-level system in Figure 2(c) can also be explained likewise, except that it 
has two additional possible radii corresponding to the same mass due to the presence of 
the second Landau level. In the one-level system in Figure 2(a), multiple values of radius 
are not observed because the magnetic field {Bo = 199.5) is such that the electrons can 
occupy only the ground Landau level. 

Finally we come to Figure 3. For each of the Figures (a), (b) and (c) the value of Epmax 
increases from the top to the bottom curve and we see that the overall radius decreases 
from the top to the bottom curve. We also note that an increase in Epmax corresponds to 
an increase in the magnetic field for a fixed Vmi which means that as the magnetic field 
strength increases the degenerate stars become more and more compact in size. The same 
trend has been observed in the case of neutron stars with a high magnetic field [26] , [27] , 
[25] (see discussion in §4.2.) Interestingly, as seen from Figure 3(a), for Ep^ax = 'iOOnieC^ 
the maximum mass of the star is even higher than that of Ep^ax = '^Orriec'^, which is 
~ 2.6Mq. Also to be noted is the fact that for a given Um, the curve corresponding to 
Epmax = 20 mgC^ covers almost the same range in mass as that of the Epmax = 200 nieC^ 
curve, while the curve for Ep^ax = '^m^c'^ covers a considerably smaller range in mass. 
The reason behind this saturation at higher Ep^ax niay be due to the fact that a Fermi 
energy ~ 2 mgC^ corresponds to very low density such that the electrons are at the most 
only mildly relativistic. By the time Fermi energy reaches a value of about 20 mgC^ (which 
corresponds to high density), the electrons have become highly relativistic, giving rise to 
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Figure 3: Comparison of the mass-radius relations, for (a) one-level system, (b) two- 
level system, (c) three-level system, when the solid, dotted and dashed lines represent 
Epmax = 2 mgC^, 20 mgC^ and 200 mgC^ respectively. In each of the three panels the j/-axis 
is in log scale. See Table 1 for details. 
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denser, more massive stars and hence further increase of Ep^ax could not bring any new 
effect in the system. 



4 Discussions 

4.1 Timescale of decay of magnetic field 

Generally the magnetic fields inside an electron degenerate star undergo Ohmic decay. 
The timescale for this is given by tohm = '^'^'^2^ , where ap is the electrical conductivity 
and L the length-scale over which the magnetic field changes. Theoretical calculations 
of Ohmic decay in isolated, cooling white dwarfs, which are generally known in nature, 
show that the magnetic field changes little over their lifetime [29], |30], [31]. Gumming 
[32] estimated a lowest order decay time as 



''ohm 



\ 1/3 / -D \ 1/2 /in ^ 



where < p > is the mean density of the star. For the three stars represented in Figure 5(d), 
above timescale turns out to be > 10^ years. Thus we can say that the magnetic fields 
of these stars do not decay significantly via Ohmic dissipation in their lifetime. However, 
while for B < 10^^ G, Ohmic decay is the dominant phenomenon, for B ~ 10^^ — 10^'^ G 
decay is supposed to take place via Hall drift and for B > 10^^ G, it undergoes ambipolar 
diffusion [33]. Note that all the decay timescales stated above are valid only for normal 
or non-superfluid matter. Hence, while they might be unrealistic for neutron stars [31] . 
these timescales are applicable to white dwarfs since they do not have a superconducting 
core [35], [36] . 

Heyl & Kulkarni [33] examined the consequences of the magnetic field decay in mag- 
netars having surface field 10^^ G to 10^^ G by using an appropriate cooling model and 
by solving the following decay equation: 

dB / 1 1 1 \ 

— = -Bi + + . (36) 

tHall I 

The strongly magnetized white dwarfs considered in the present work have central mag- 
netic field strengths ~ 10^^ G, which is comparable to the surface field strengths of 
magnetars. It is then likely that the magnetic fields in both these cases would undergo 
similar decay mechanisms. The decay of such strong fields is dominated by ambipolar 
diffusion [33]. However, for a typical initial field strength of 10^^ G, the magnetic field 
remains nearly constant up to about 10^ years and in the next 100 years its value de- 
creases by at most an order of magnitude [23j. As discussed in §4.3 below, it is the central 
magnetic field which is crucial for the super-Ghandrasekhar mass of the white dwarfs. 
Thus applying the above results to these magnetized white dwarfs, we can conclude that 
the central field will not decay appreciably for a long period of time. 

An alternate scenario could arise if the magnetized white dwarfs are accreting. In 
this case, the heat generated due to accretion decreases the electrical conductivity of the 
surface of the star, causing a faster decay of the (surface) magnetic field due to a reduced 
Ohmic decay timescale. However, as the mass of the star increases, it becomes more 
compact and the current carrying accreted material is pushed deeper into the star. Since 
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conductivity is a steeply increasing function of density, the higher conductivity of the 
denser inner region of the star will again slow down further decay of the magnetic field 

Since we are working with magnetic field strengths B > Be, one might be concerned 
about the process of electron-positron pair creation (via Schwinger process) at the expense 
of magnetic energy, which might lead to a reduction of the field strength. However, Canuto 
and Chiu [3H], |39], |10] and also Daugherty et al. [H] showed that it is impossible to 
have spontaneous pair creation in a magnetic field alone, irrespective of its strength. Now, 
from Maxwell's equations in a steady state, we have 

V X B = — j (37) 
c 

and Ohm's law states 

E = — , (38) 

when E is the electric field and j the current density. Thus for an electric field to be 
generated, the magnetic field must vary with space as seen from equation ( 137|) . However, 
in this work we have chosen a constant magnetic field, which leads to j = and hence 
E = from equation ( l38ll . Now the main effect of the magnetic field, which is to give 
rise to a mass exceeding the Chandrasekhar limit, is restricted to the high density region, 
where the field remains essentially constant. Thus our choice of a constant magnetic field 
is justified (see §4.3 for detailed discussion). Moreover since the magnetic field in the 
degenerate star is generated due to the flux freezing phenomenon, it incorporates the 
fact that the conductivity aE is very large. Thus even if the magnetic field is highly 
inhomogeneous, from equation fl55]) we see that the electric field generated would again 
be negligible. Thus the electron degenerate stars in this work are magnetically dominated 
systems, i.e., both the magnetic field strength is very high {B > Be) and the electric field 
strength is negligible. Recently, Jones |12] calculated the cross-section for photon-induced 
pair creation in very high magnetic fields and has arrived at the result that, there is a 
rapid decrease of the pair creation cross-section at B > Bf.. Hence one can ignore the 
effect of pair creation in reducing the magnetic field strength in these stars. 

4.2 Comparison with Chandrasekhar 's results 

In Figure 4 we put together our results along with Chandrasekhar 's result for non- 
magnetic white dwarfs to obtain a more complete picture. Figure 4(a) shows the equation 
of state obtained by Chandrasekhar and those corresponding to the one (z/^ = 1), two 
{um = 2) and three = 3) level systems (same as in Figure 1(b)). Interestingly, the 
equation of state for z/^ = 20 almost grazes Chandrasekhar 's equation of state, except 
the appearance of a series of kinks. This clearly shows that as the magnetic field strength 
decreases, or equivalently as the maximum number of occupied Landau levels increases, 
the equation of state approaches Chandrasekhar 's equation of state. 

Figure 4(b) represents the mass-radius relations corresponding to the equations of 
state shown in Figure 4(a). The mass-radius relation for z/^ = 500 (500 Landau levels) is 
also shown. We observe that as the magnetic field strength decreases or as the maximum 
number of occupied Landau levels increases, the mass-radius relation approaches the non- 
magnetic relation and one recovers Chandrasekhar 's mass limit. 
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Figure 4: Comparison with Chandrasekhar's non-magnetic results for Epmax = 20 mgC^. 
(a) Equations of state - the sohd hne represents Chandrasekhar's equation of state. The 
dot-dashed, dotted and dashed hnes represent the one-level = 1), two-level — 2) 
and three-level — 3) systems respectively. The equation of state for Um — 20 is also 
shown, which appears as a series of kinks on top of the solid line, (b) Mass-radius relations 
- the vertical line marks the 1.44Mq limit and the solid line represents Chandrasekhar's 
mass-radius relation. From top to bottom the other lines represent the cases for i/^ = 
500 , 20 , 3 , 2 and 1 respectively (the |/-axis is in log scale), (c) Density as a function of 
radius inside a non-magnetized star (solid line) and a star having — 1 (dot-dashed 
hne), both having the same central density pc = 5 in units of 2 x 10^ gm/cc. 
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Figure 4(b) also shows that the stars of a given mass become more and more compact 
in size as the magnetic field strength increases. In order to explain this behavior we 
resort to Figures 4(a) and (c). Let us for simplicity look at Chandrasekhar's equation of 
state (solid line) and that of the one-level system (dot-dashed line), which corresponds 
to B = 199. 5Bc, in Figure 4(a). We notice that at low densities, the dot-dashed line lies 
below the solid line and at higher densities, > 2, the dot-dashed line lies above the 
solid line. In other words, the equation of state for the one-level system is softer than 
Chandrasekhar's equation of state at low densities, which means that the pressure does 
not rise with density as rapidly as that in Chandrasekhar's case. This trend reverses 
at higher densities and the equation of state for the one-level case becomes stiffer than 
that of Chandrasekhar's. Now, matter with a softer equation of state is less efficient in 
counteracting gravity and hence stars made out of such matter will be more compact in 
size. Keeping this in mind we now look at Figure 4(c). Figure 4(c) shows the variation 
of density with radius, within a non-magnetized star and a magnetized star with 1/^ = 1 
(one- level), both having the same central density {pc = 5). We mention here that for 
ease of explanation we have chosen a high central density such that the equations of state 
for both the stars cover almost the entire range of density. We observe that for a given 
radius, density for the magnetized star is higher than that of the non-magnetized star 
for a large range, from = Pc = 5 to p/) ~ 0.3. But at very low densities, pd < 0.3, 
the density of the magnetized star sharply falls to zero due to smaller pressure, leading 
to a smaller star {R = 6.4 x 10'' cm) than in the non- magnetized case {R = 1.3 x 10^ 
cm). From the analysis of the equations of state for these two stars, we can say that if 
the equation of state of the magnetized star would have remained stiffer than the non- 
magnetized star throughout, then the magnetized star would have a larger radius, as in 
this case the pressure would be more efficient in counteracting the gravitational collapse 
of the star. But this is not true. At very low densities the equation of state suddenly 
becomes softer for the magnetized star and the pressure is not able to counteract gravity 
efficiently, causing the star to collapse rapidly. Hence the density goes to zero very rapidly, 
causing the star to have a (much) smaller radius. 

We also note from Figure 4(b) that as the magnetic field strength increases, the prob- 
ability that the stars will have masses exceeding the Chandrasekhar limit increases. From 
Figure 4(c) we see that the density of the magnetized star is much higher than that of 
the non-magnetized star, except for a very small range of low densities. Thus calculating 
the total mass of the stars from equation (12^ . we obtain a much higher value for the 
magnetized star. 

4.3 Choice of constant magnetic field 

Figures 5(a), (b) and (c) show the variation of mass as a function of density within 
a magnetized electron degenerate star for three different magnetic field strengths. In all 
the cases we note that, by the time the density falls to about half the value of the central 
density, the mass has increased significantly, crossing the Chandrasekhar limit (indicated 
by the horizontal line) soon after. Hence, although we have considered a constant magnetic 
field, their effect is restricted to the high density regime, where the field remains essentially 
constant in reality. Hence we can also interpret this constant magnetic field as the central 
magnetic field of the star. This would be more clear from the description of the variation 
of magnetic field given by [26], [13], which show that an inhomogeneous magnetic profile 
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Figure 5: Mass as a function of density within an electron degenerate star with Bjj 
(a) 199.5 (dashed hne), (b) 99.75 (dotted hne), (c) 66.5 (dot-dashed hne), for Epmax — 
20 meC^. (d) Variation of mass as a function of radius within the stars in (a), (b) and (c). 
The horizontal line indicates the 1A4Mq limit. The radius is in units of 10^ cm. See 
Table 1 for details. 
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in a compact star could be such that the magnetic field is nearly constant throughout 
most of the star and then gradually falls off close to the surface (see Figure 5(b) in [26j). 
Thus choosing an inhomogeneous magnetic profile would not affect our main finding that 
the Chandrasekhar mass limit can be exceeded for high magnetic field strengths. 

4.4 Choice of non-relativistic equation of hydrostatic equiUb- 
rium 

Figure 5(d) shows the variation of mass as a function of radius within the stars rep- 
resented in Figures 5(a), (b) and (c). All these stars have a total mass ~ 2.3Mq and 
hence their Schwarzschild radius Rg = = 6.8 x 10^ cm. Now, general relativistic 
effects usually start becoming important at a radius < lORg, i.e. ^ 0.068 for the 
above mentioned stars. However, from Figure 5(d) we see that for all the three stars 
the contribution to the mass from a radius < lORg, i.e. the central region, is negligible. 
Contribution to the mass rather effectively starts from a radius vd > 0.13 (~ 20Rg) and 
the Chandrasekhar limit is crossed at Rg ~ 66, 95 and 118, for the stars represented by 
the dashed, dotted and dot-dashed lines, respectively. Thus significant contribution to the 
total mass of the stars comes from a region well beyond the regime of general relativity. 
Hence one need not necessarily consider the Tolman-Oppenheimer-Volkoff equation and 
our choice of the non-relativistic equation of equilibrium, equation (JT9l) . is justified. 

4.5 Unstable branch of the mass-radius relations 

Figure 6(a) shows the variation of density as a function of radius inside the degenerate 
stars having Epmax = 20 mgC^ and B = 99.75i?c- If we look at the equation of state 
for the two-level system with Epmax = 20 mgC^, i.e. the dotted line in Figure 1(b), we 
see that both pu = 1.5 and 2 lie below the density at which transition takes place from 
ground Landau level to first Landau level. For stars with these central densities, the 
pressure rises monotonically with density throughout and one forms a stable star. From 
the corresponding curves (solid and dashed) in Figure 6(b) also we note that the sound 
speed (cg) varies smoothly with radius, reaching a maximum value at the center. Both 
the mass and radius of these stars increase with central density and hence they lie on the 
uppermost branch of the mass-radius relation in Figure 2(b). 

Now pd = 2.5 lies on the plateau following the kink in the equation of state. A star 
having this central density will have within it a zone where pressure does not steadily 
increase with density, but is nearly constant. Thus such a star will tend to collapse faster 
under gravity and hence will have a smaller radius. From the dotted line in Figure 6(a) 
we see that the density of such a star falls steeply to zero compared to the pu = 2 case 
indicated by the dashed line. The corresponding curve in Figure 6(b) shows a peak in Cg 
at the radius where pressure starts becoming independent of density. Since mass of the 
star is calculated using equation (129]) . they have smaller masses too due to smaller size. 
These stars lie on the middle branch of the corresponding mass-radius relation. Since 
these stars consist of an unstable zone mentioned above, we argue that they constitute 
the somewhat unstable branch in the mass-radius relation. 

Finally, let us choose a larger density, say pd = 7, which lies on that portion of the 
equation of state above the kink, where pressure again rises monotonically with density, 
but not as steeply as that in the regions around po = 1-5 and 2. Although these stars 
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Figure 6: (a) Density and (b) sound speed in units of c, as a function of radius, inside 
the stars having Epmax = 20 mgC^ and B = 99.75i?c- The sohd, dashed, dotted and dot- 
dashed hues represent the stars with central densities 1.5, 2, 2.5 and 7, in units of 2 x 10^ 
gm/cc, respectively. 



also have an unstable zone inside them, but softer pressure in a large range of density, 
which together cause them to have a smaller radius, their central density is very high. As 
a result their mass starts increasing again, as can be inferred from the dot-dashed line 
in Figure 6(a). The again shows a sharp peak indicating the appearance of unstable 
zone within the star. These stars constitute the bottommost branch of the mass-radius 
relation. 

4.6 Neglecting Coulomb interactions 

The distance between nuclei in highly magnetized electron degenerate stars could be as 

2 

small as 50 Fermi. The Coulomb repulsion energy between electrons — , at the separation 
of 50 Fermi, is of the order of 4 x 10~^ ergs which is quite less than the rest mass energy 
of an electron mgC^ ~ 8 x 10"'' ergs. There can be Coulomb interaction between electrons 
and the ions, which is given by Commonly electron degenerate stars consist of 

helium, carbon, oxygen, etcetera, so that Z can have a value of 10 at the most. Hence, 
the Coulomb interaction energy would still be less than, or at most the same order as, the 
rest mass energy. Thus we can neglect the effects of Coulomb interaction for the present 
purpose. 

4.7 Anisotropy in pressure due to strong magnetic field 

The strong magnetic field causes the pressure to become anisotropic [27], [44], [45] . 
The total energy momentum tensor due to both matter and magnetic field is to be given 

by 

where, 

= e^u^u'' - PUg"" - W-un (40) 
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when, em = Ce is the matter energy density given by equation (fT5l) and Pm = Pe is the 
matter pressure given by equation (IT7I) . The first term in equation ( HTj) is equivalent to 
magnetic pressure, while the second term gives rise to the magnetic tension. If B is along 
the 2;-axis, then we have 
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Thus we see that pressure becomes anisotropic. The total pressure in the perpendicular 
direction to the magnetic field is given by 
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and that in the parallel direction to the magnetic field is given by 
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Now the parallel pressure becomes negative if the magnetic pressure exceeds the fluid 



pressure. In order to understand this effect, we write the component of T 
z-axis as 
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The second term —B'^/An corresponds to an excess negative pressure or tension along the 
direction to the magnetic field. Thus the total energy momentum tensor can be written 

as 



rp/ll/ 



. +M1 











p +^ 









p + 

-'■ml 





Stt 



{Pm + 







Bl 

Stt . 



Ml 

Atx 



(46) 



The strong magnetic field also reveals anisotropy due to magnetization pressure 
Hence, actually the pressure in the perpendicular direction to the magnetic field is given 
by 

52 



Pm + — -MB, 



where M. is the magnetization of the system, which is given by 
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(47) 
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However, for the magnetic fields considered in the present work exhibiting super-Chandrasekhar 
masses, i?2/87r >> M.B. Therefore, we do not include magnetization term in the pressure 
which does not affect the result practically for the present purpose. 
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Here we refer to the work by Bocquet et al. [16] , which models rotating neutron stars 
with magnetic fields, by using an extension of the electromagnetic code used by Bonazzola 
et al. ^47j . They observed that the component of the total energy momentum tensor along 
the symmetry axis becomes negative (equivalent to T^^ < in our case), since the fiuid 
pressure decreases more rapidly than the magnetic pressure away from the center of the 
star. This happens because the combined fiuid-magnetic medium develops a tension. As 
a result of this magnetic tension, the star displays a pinch across the symmetry axis and 
assumes a fiattened shape. A similar effect is expected to occur in our work, where the 
magnetic tension will be responsible for deforming the magnetized white dwarf along the 
direction to the magnetic field and turns it into a kind of oblate spheroid. Hence, one 
should be cautious before considering equation (fT9!) which is applicable for a spherical 
star |I5] . 

However, in this work we consider a constant magnetic field, since our interest is to see 
the effect of the magnetic field of the central region of the white dwarf, where the field is 
supposed to be (almost) constant (see §4.3). Thus even if we use either the parallel or the 
perpendicular pressure in the hydrostatic equilibrium equation f|T9|) . B does not appear 
explicitly in the equation [dB/dr = 0) — only the gravitational field will be modified due 
to deformation. Hence, it is still possible to have super-Chandrasekhar mass white dwarfs 
- only they will be deformed in shape due to the strong field, which might even render a 
more massive white dwarf (as is discussed in Appendix A). 

5 Summary and Conclusions 

We have studied the effect of high magnetic field on the equation of state of purely 
electron degenerate matter at zero temperature. In the equation of state, we have consid- 
ered only the electron degeneracy pressure modified by the strong magnetic field. We have 
focused on those Landau quantized systems in which the maximum number of Landau 
level(s) occupied is/are one, two or three, which we have named to be one-level, two-level 
and three-level system respectively. 

We have found that whenever a lower Landau level is completely filled and the next 
higher level is to be filled, a kink appears in the equation of state, followed by a plateau - 
which is a small region where the pressure becomes nearly independent of the density. The 
one-level system, which has only the ground Landau level filled, has no kink, the two-level 
system has one kink at the ground to first level transition and the three-level system has 
two kinks, one at the ground to first and the other at the first to second-level transition. 
We have studied each of these systems at three maximum Fermi energies Epmax = 2 Tn^c^, 
20 nieC^ and 200 mgC^ and obtained the mass-radius relations of the corresponding stars. 

The mass-radius relations show turning point (s), denoted by the decrease of both mass 
and radius, which correspond(s) to the kink(s) in the equation of state. The most inter- 
esting result obtained is that there are possible stars found on the mass-radius relations 
whose mass exceeds the Chandrasekhar limit. They could be potential magnetized white 
dwarfs. The maximum mass obtained is about 2.3 — 2.6M0 and is seen to occur for various 
combinations of central density, magnetic field strength and maximum number of occu- 
pied Landau levels. Interestingly such super-massive white dwarfs have been suggested 
to be the most likely progenitors of recently observed Type la supernovae [20], [21]. Out 
of the equations of state considered in the present work, the system with the lowest mag- 
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netic field which gives rise to this mass is the three-level system with Epmax = 20 mgC^, 
Bp, = 66.5 (or B = 2.94 x 10^^ G) and the corresponding central density being 2.2 x 10^ 
gm/cc. 

The nature of the mass-radius relations is governed by the fact whether the system 
is one- level, two-level or three- level and is independent of the value of Epmax- However, 
Epmax determines how relativistic the system is. For instance, the Chandrasekhar mass 
limit is not exceeded for a low Epmax (say = 2meC^), no matter what the central density 
is. We have however observed that as Epmax increases, which corresponds to an increase 
in the magnetic field strength, the degenerate stars become more compact in size. 

As discussed, the minimum magnetic field required to have a 2.3M0 degenerate star 
is -B = 2.94 X 10^^ G. The magnetic field of the original star of radius Rq, which collapses 
into the above degenerate star of radius ~ 10^ cm, turns out to be ~ 6 x 10^ G, based 
on the flux freezing theorem. Existence of such stars is not ruled out [15]. However, the 
anisotropy in pressure due to strong magnetic field causes a deformation in the white 
dwarfs which adopt a flattened shape. This effect of flattening leads to more massive 
white dwarfs, even at relatively lower magnetic field strengths. 

One might wonder as to why highly magnetized white dwarfs have not yet been ob- 
served. A plausible reason could be that the surface magnetic field is being screened due 
to some physical processes. For instance, if the white dwarf is in a binary system and is 
accreting matter from its companion, as is proposed for Type la supernovae progenitors, 
then the plasma that is being deposited on the surface of the star could induce an opposite 
magnetic moment. This would result in a reduction of the surface field strength. How- 
ever, the central magnetic field strength, which is presumably unaffected by the above 
processes, could be several orders of magnitude higher than the surface field. Indeed as 
seen in §4.3, it is the central field which is crucial for exceeding the Chandrasekhar mass 
limit. 
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A Appendix 

In order to estimate the effect of deviation from spherical symmetry due to the mag- 
netic field we have performed a few calculations. If the magnetic field is very strong, then 
the magnetic tension will flatten the star along the direction to the field (as discussed in 
§4.7). If we consider the white dwarf to be an oblate spheroid with the z-axis being the 
symmetry axis, then its equation is given by 

^ + ^ = 1, (49) 

where, = + y^, the semi- axis a is the equatorial radius of the spheroid and c is the 
distance from center to the pole along the symmetry axis (c < a for an oblate spheroid). 
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The equatorial force balance equation can be written as 
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(50) 



where, = r^^ + . The above equation for hydrostatic equilibrium has to be supple- 
mented by an equation to determine the mass of the star. Now, 

dM = pdV (51) 

and the volume element for an oblate spheroid is given by 

dV = Tireq{zfdz. (52) 

For a simpler visualization, one can also assume a cylindrical geometry and the equation 
for the mass can be given by 

dM , 

- — = 2nreghp, (53) 

where h denotes an average height of the cylinder (assuming that the density does not 
change appreciably with z) and is a parameter that quantifies the degree of flattening. 
At high field strengths, the white dwarf will be more flattened and h will be small. At 
low field strengths it is likely that the star will be less flattened and h will have a larger 
value. Keeping this in mind we solved the above equations and obtained the mass-radius 
relations for some of the cases as shown in Figure [71 

The one-level and two-level systems have strong magnetic fields and the corresponding 
mass-radius relations for the spherical case are denoted by the solid lines in Figures 7(a) 
and (b). From the dashed lines in Figures 7(a) and (b) we see that the flattened white 
dwarfs could have much higher masses than the perfectly spherical ones. Interestingly we 
note from the dashed lines of Figures 7(c) and (d) that even for much lower magnetic field 
strengths {B = 4.4 x 10^^ G for z/„ = 20 and 1.8 x lO^^ G for = 500), the Chandrasekhar 
mass limit is exceeded even if one includes the corresponding reduced flattening effect - 
these white dwarfs were sub-Chandrasekhar for the spherical cases. One can also see that 
the radius of the stars represented by the dashed lines is larger than that represented by 
the solid lines. One must note here that for the dashed lines this radius is the equatorial 
radius, which will automatically be larger than the spherical radius as a consequence of 
the flattening effect. 

Hence, the effect of flattening leads to a more massive star. The effect is very similar 
to the flattening due to centrifugal force in rapidly rotating stars, which are known to 
be more massive than their slow-rotating counterparts (e.g. [46J). These are also shown 
to have larger mass in presence of high magnetic field. Thus the strong magnetic field 
is responsible for the increasing the mass of the white dwarfs, while the deformation (or 
flattening) of the white dwarf due to the field further adds on to the mass. Therefore, our 
estimate of the mass of the white dwarf, in fact just sets a lower bound. More interestingly, 
flattening effects due to magnetic field will render super-Chandrasekhar white dwarfs even 
at a smaller magnetic field - such relatively low magnetized white dwarfs are more probable 
in nature. 
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(a) (b) 
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Figure 7: Mass-radius relations with Epmax — 20 mgC^ for (a) one-level system — 1), 
(b) two-level system (z/^ = 2), (c) twenty-level system = 20) and (d) five hundred- 
level system (z/^ = 500). Here Md is the mass of the white dwarf in units of Mq and 
Reg is the equatorial radius of the white dwarf in units of 10® cm. All the solid lines 
represent the mass-radius relations for the cases, if the stars would have been spherical. 
The dashed hues in (a) and (b) represent the mass-radius relations for highly flattened 
(strongly magnetized) white dwarfs, while the dashed lines in (c) and (d) represent the 
mass-radius relations for less flattened (relatively weakly magnetized) white dwarfs. In 
all the four panels the y-axis is in log scale. 
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